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Abstract 

Molecular techniques have received increased interest in many ecological studies. The 
ability of such techniques to determine population structure, gender, and even individuals has 
provided valuable insight to effective management of wildlife populations. A major drawback 
of such applications in wildlife studies is the invasiveness of traditional techniques of DNA 
collection, many of which require capture of individuals. In oviparous chordates, the use of 
eggshells as sources of genetic material would mitigate the impact of DNA collection on the 
species of interest. Here we present procedures and protocols used to extract genomic DNA 
from Rio Grande wild turkey (Meleagris gallopavo intermedia) eggs collected at two study sites. 
While DNA concentrations were highly variable (37.27 ng/pl ±5.12 SE), the below procedures 
resulted in reliable genetic data for both gender identification as well as microsatellite analysis 
using previously developed primers. Based on our findings, eggshell remains appear to have 
enormous potential as a noninvasive source of DNA. 

Key words: DNA, eggshell, Meleagris gallopavo , noninvasive, PCR, sex ratio, wild 
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Introduction 


Genetic analyses of wildlife populations have 
contributed greatly to our understanding of popula¬ 
tion structure and dispersal, conservation biology, 
and wildlife forensics (Haig 1998; Parker et al. 1998). 
However, the most commonly used sources of wildlife 
genetic material, blood or tissue samples, are invasive, 


requiring the capture and restraint of individual animals 
(Piggot and Taylor 2003). Invasive techniques may 
negatively affect the species of interest either physically 
or behaviorally and are often logistically problematic 
(Spraker et al. 1987; Nicholson et al. 2000). Such 
concerns, coupled with increasing technologies allow- 
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ing DNA extraction from dilute or small samples, have 
lead to exploration of noninvasive alternatives for DNA 
collection (Taberlet et al. 1999). 

Noninvasive genetic sampling techniques have 
played larger roles in wildlife ecology in recent years 
(Taberlet etal. 1999; Pearse et al. 2001). DNAextract- 
ed from hairs has been used to estimate population size 
and gene diversity in black bears, Ursus americanus 
(Boersen et al. 2003; Triant et al. 2004). DNA extracted 
from feces has been used to determine genera in canids 
(Reed et al. 2005) and felids (Farrell et al. 2000) and for 
gender determination in several species (e.g., Reed et 
al. 1997; Ortega et al. 2004). While both hair and feces 
have been shown effective as noninvasive techniques 
for mammalian species, noninvasive techniques in the 
majority of oviparous species are less prevalent. 


Materials 

Field Collection of Eggshell and Blood 
Samples - Rio Grande wild turkey eggs were collected 
at two study sites in Collingsworth, Cottle, and Donley 
counties, Texas, (Texas Panhandle) in conjunction 
with ongoing Rio Grande wild turkey research at 
Texas Tech University. Only egg remains that could 
be classified as unique (i.e., collect only blunt or 
pointed ends at a given nest site) were collected. Each 
clutch was assigned a unique number, and each egg 
from that clutch was numbered and placed in separate 
plastic storage bags. All samples were stored at 4°C 
within 12 hours of collection and transported from 
field sites to the laboratory in coolers. For verification 
purposes, blood samples (2-3 drops) were collected 
from the brachial vein of adult and yearling wild 
turkey during capture and placed directly in lysis 
buffer (Longmire et al. 1997). 

Storage of Eggshell Matena/.-Chorioallantoic 
membrane (-1.5 cm 2 ) from each egg sample was 
excised and placed in a sterile petri dish with 1 ml of 
lysis buffer (Longmire et al. 1997). To facilitate the 
appropriate lysing and subsequent digestion steps, the 
membrane was ground into smaller pieces using a ster¬ 
ile single-edge razor blade. The ground membrane/ 
lysis mixture was placed into a 5 ml polypropylene 
tube and brought to a total volume of 2 ml with fresh 


In avian species, feathers have traditionally been 
used for DNA extraction (Taberlet and Bouvet 1991), 
with recent studies acquiring resolutions to address 
population size, structure, and parentage (McCracken 
et al. 2001; Segelbacher and Storch 2002; Rudnick et 
al. 2005). Studies also suggest eggshell remains may 
be a valuable source of genetic material with distinct 
advantages for wildlife managers (Pearce et al. 1997; 
Strausberger and Ashley 2001). Using established 
primers (D’Costa and Petitte 1998; Latch et al. 2002) 
and sampling protocols (Longmire et al. 1997) based 
on more traditional collection methods (e.g. tissue, 
blood), we attempted to determine the efficiency of 
DNA extraction from field collected Rio Grande wild 
turkey ( Meleagris gallopavo intermedia) eggs for 
gender determination, as well as its logistic feasibility 
for wildlife studies. 


and Methods 

lysis buffer. The work area and all instruments were 
first washed with 70% EtOH and then wiped clean 
with sodium hypochlorite (6% bleach) to reduce the 
risk of cross contamination. 

DNA Isolation-D NA extraction was carried 
out using the method described by Longmire et al. 
(1997) with modifications. Extractions began by 
transferring approximately half of the original mem¬ 
brane/lysis solution from the 5 ml tube into a fresh 
15 ml Falcon tube. The volumes in both tubes were 
brought back up to 1.5 ml with fresh lysis buffer. The 
remaining original sample was archived in the Natural 
Science Research Laboratory (NSRL) ( http://www. 
nsrl.ttu.edu/) at the Museum of Texas Tech Univer¬ 
sity. To the Falcon tube, 30 jul of Proteinase K, at a 
concentration of 10 mg/ml, were added to the 2 ml 
membrane/lysis solution. Digestion was carried out 
overnight (-15 hrs) in a 37°C incubator with slow 
rotation. After overnight digestion, the liquid phase 
was decanted into a new 15 ml tube, and the remain¬ 
ing solid material was discarded. An equal volume 
(-1.5 ml) of TE buffer saturated phenol was added to 
the digested membrane solution. The phenol wash 
was carried out in a 37°C incubator with slow rota¬ 
tion for approximately 30 minutes. The aqueous and 
organic phases were separated by centrifugation at 
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2000 rpm for 10 minutes. The aqueous phase (top) 
was removed and placed into a clean 15 ml tube, and a 
second phenol wash was conducted. After the second 
phenol wash, a final extraction wash using Phenol: 
Chloroform:Isoamyl alcohol at a mixture of 25:24:1 
as outlined above was used. 

The remaining aqueous phase was put into dialy¬ 
sis tubing (MWCO 12.14,000) strips approximately 
7 cm in length. Dialysis was carried out in a IX TE 
buffer solution at 4°C for 72 hours, changing the buffer 
every 3-4 hours. Total genomic DNA (gDNA) was 
stored at 4°C for immediate use or frozen (-20°C) for 
archival purposes. The quantity of isolated DNA was 
measured using UV spectrophotometric analysis of 
A260/A280 ratio on an ND-1000 spectrophotometer 
(NanoDrop Technologies). The isolated DNA was 
visualized using agarose gel electrophoresis (0.8%) 
and stained with ethidium bromide to assess quality 
(Fig. 1). 

Sex Identification using PCR -The PCR reaction 
required a substantial amount of template due to the 
small quantities of DNA isolated from the eggshells. 
The PCR protocol and thermal profile was modified 
from D’Costa and Petitte (1998) which was developed 
for sex identification of turkey embryos. Final con¬ 
centrations of the optimized multiplex-PCR (M-PCR) 
conditions were: 0.2 mM of each dNTPs, 0.5 pM of 
each primer (ATP synthase gene and PstI primer sets), 
0.8 units of Taq DNA polymerase (Promega), 2 mM 
MgCl 2 (Promega) in a IX Mg free buffer (Promega) 
with 0.25 mg/ml of bovine serum albumen (bsa) 
and 150 ng of isolated DNA. All reactions were 


performed at a final volume of 10 pi in 0.2 ml thin 
wall PCR tubes. The amplifications were performed 
using a PTC-200 Gradient Cycler (MJ Research) 
with the following thermal profile: initial denature at 
94°C for 2 minutes; 40 cycles of 94°C for 30 seconds 
(denature), 51°C for 30 seconds (annealing) and 72°C 
for 40 seconds (extension); final extension of 72°C 
for 10 minutes. The isolated gene products for both 
the ATP synthase (250 bp) and PstI repeat (177 bp) 
were visualized using agarose gel electrophoresis 
(2%) and stained with ethidium bromide to visualize 
the amplicons by band size discrepancy under a UV 
light source. 

Microsatellite Amplification - We amplified and 
sequenced the WT75 and WT38-2 loci from 4 indi¬ 
viduals (Latch et al. 2002). We used 2 samples (one 
male and one female) of DNA isolated from whole 
turkey blood collected during winter capture events 
(2005) and 2 samples (one male and one female) of 
DNA isolated from turkey egg shell remains. PCR 
conditions followed those from Latch et al. (2002), 
and the final concentrations for PCR conditions were: 
0.2 mM of each dNTPs, 0.5 pM of each primer, 0.65 
units of Taq DNA polymerase (Promega), 1.5 mM 
MgCf (Promega) in a IX Mg free buffer (Promega) 
with 10 ng of isolated DNA. All reactions were per¬ 
formed at a final volume of 50 pi in 0.2 ml thin wall 
PCR tubes. The amplifications were performed using 
a PTC-200 Gradient Cycler (MJ Research) with the 
following thermal profile: initial denature at 95°C for 2 
minutes; 30 cycles of 95°C for 30 seconds (denature), 
55°C (for WT75) or 59°C (for WT38-2) for 30 seconds 
(annealing), and 72°C for 40 seconds (extension); final 
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Figure 1. Total genomic DNA isolated from the chorioallan¬ 
toic membrane of Meleagris gallopavo intermedia. Lanes 
1-9 represent nine individual samples from one clutch. Lane 
M is a 1 kilobase (kb) ladder. 
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extension of 72°C for 10 minutes. The isolated gene 
products for both the WT75 locus (278-298 bp) and 
WT38-2 locus (103-119 bp) were visualized using 
agarose gel electrophoresis (2%) and stained with 
ethidium bromide to visualize the units by band size 
discrepancy under a UV light source. 

All PCR products were purified before sequenc¬ 
ing reactions were performed using the QIAquick* 
Spin Kit (QIAGEN) according to manufacturer’s 


protocol. Both strands of the individual amplicons 
were sequenced (Sanger et al. 1977) on an ABI 3100 
Genetic Analyzer (ABI) according to ABI protocol. 
Sequences were proofed, assembled and aligned using 
ContigExpress and AlignX (respectively) in Vector NTI 
Suite 7 (InforMax, Invitrogen life science software). 
Alignments were made based on sequence data depos¬ 
ited in GenBank (NCBI) for the WT75 locus (accession 
number AF434907) and for the WT38-2 locus (acces¬ 
sion number U79365) (Latch et al. 2002). 


Results 


Sex Identification using PCR .-From eggs col¬ 
lected no more than 9 days post-hatch, we assigned 
sex identifications for 63 individuals from a total of 
7 unique clutches over two separate study sites (Ap¬ 
pendix 1). Of these 63 individuals, 35 were identified 
as males, showing a single band, and 28 were females, 
showing two bands (Fig. 2). The Matador site (n = 20) 
consisted of 70.0% male individuals, and the Salt Fork 
site (n = 43) 48.9% males. To assess accuracy of the 
technique, we correctly identified 6 individuals (n = 3 
males and 3 females) using DNA isolated from blood 
of known male and female turkeys (Fig. 3). 


Microsatellite Amplification .-We were able to 
successfully amplify both the WT75 (Fig. 4) and the 
WT38-2 loci from both eggshell and blood DNA. The 
repetitive elements for the Salt Fork site individuals 
were [TG] 4 and [TCG] 7 for TK136464 and [TG] n for 
TK136466. The repetitive elements for the Matador 
site individuals were [GT] . and [CT] 10 for M307-2 and 
[GT] n and [CT] 7 for M307-4. For the WT38-2, the 
repetitive element for both study sites was [CA] |5 . 


1234567 89M 



Figure 2. Sex identification of turkey egg shell remains by 
multiplex PCR (M-PCR) using total genomic DNA. Lane M 
is a 1 kb ladder. All individuals amplify the 250 bp turkey 
specific ATP synthase gene, and females are determined by 
the 177 bp PstI repeat. Male samples are shown with one 
band in lanes 2-3 and 6-8. Female samples are shown with 
2 bands in lanes 1 and 4-5. Lane 9 is a no template control 
reaction. 
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Figure 3. Sex identification by multiplex PCR (M-PCR) using 
total genomic DNA isolated from blood of adult birds of known 
sex. Lane M is a 1 kb ladder. Female samples are shown with 
two bands in lanes 1 -3. Male samples are shown with one band 
in lanes 4-6. Lane 7 is a no template control reaction. 
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Figure 4. Amplification and PCR purification of the WT75 
locus using total genomic DNA isolated from blood and 
egg shell remains. Lane M is a low mass ladder. DNA 
from blood is in lanes 1-2. DNA from egg shell remains 
is in lanes 3-4. 


Discussion 


Eggshell remains are a logistically feasible, non- 
invasive source of genetic data. While our study did not 
empirically address DNA degradation over time, viable 
DNA was obtainable over a rather large temporal win¬ 
dow (9 days) between nest hatch and collection of eggs. 
Environmental conditions such as heat, sunlight, rain, 
or insects may have initially affected DNA quantity or 
quality. Post collection, time to suspension appeared 


to influence DNA quantity as evidenced by our trend 
line (Fig. 5). In radio telemetry studies, eggs could 
be collected with little additional effort and provide a 
wealth of genetic information on a range of questions, 
possibly even providing genetic data from the mother 
(Strausberger and Ashley 2001). DNA from eggshells 
may prove to be an easily collected and stored ware¬ 
house of genetic material. 
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Figure 5. A trend line illustrating degradation of final product DNA 
concentrations between collection time and suspension in a lysis buffer 
solution (i? 2 =0.48, n=6 3). 


For wildlife managers, sex ratios at hatch are a 
mystery in many game bird species. Several studies 
suggest equal sex ratios may be inaccurate in some 
avian species (Hamilton 1967; Rosenfield et al. 1996; 
Sheldon 1998), yet wild turkey population models 
typically assume a 50/50 sex ratio (Healy and Powell 
1999). While Healy and Powell (1999) suggest skewed 
sex ratios may influence models, the data presented here 
are the only datasets known to the authors regarding 
wild turkey sex ratios at birth. In fact, the paucity of 
sex ratio data on many game species is indicative of the 
difficulty associated with capturing cryptically-colored, 
sexually monomorphic, precocial young. Further, 
direct capture of young may invite sexually-biased 
capture probabilities and influence other population 
estimates associated with study goals (Spraker et al. 
1987; Colwell etal. 1988). However, eggshells may be 
easily collected, inhibit repetitive sampling, and main¬ 
tain unbiased population estimates. While the authors 
foresee many potential uses for eggshell DNA (e.g., 
dump-nesting documentation, parentage), none may 
be as beneficial as the ability of eggshells to produce 
estimates of secondary sex ratios. 


We believe that we have a straightforward and 
accurate method for wildlife and molecular biologists 
to use wild turkey egg shell remains as a source of 
genetic material. Our primary goal was to test whether 
or not sex identifications could be made using DNA 
isolated from egg shell remains (i.e., chorioallantoic 
membrane). The use of a turkey specific ATP synthase 
gene may allow for accurate gender determination 
even in possibly contaminated samples (e.g., avian 
predation). Amplifying the ATP synthase gene in a 
multiplex reaction with the female specific PstI element 
greatly reduces the cost of this method. An important 
note when using this protocol is the amount of DNA 
necessary to run the sexing PCR. The range of DNA 
concentrations was 0.6 to 214.3 ng/pl (x = 37.27 ± 
5.12) for samples used in this study. This indicates 
that some reactions, limited by the reaction volume, do 
not actually have 150 ng of DNA template. While we 
had a number of reactions work using approximately 
100 ng of starting material, 150 ng of starting material 
is optimal for this protocol and should be the target 
concentration. 
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We wanted to further assess the quality and 
capabilities of DNA from egg shell remains for other 
molecular techniques. A number of studies on the ge¬ 
netic variation in wild and domestic turkeys using short 
tandem repeats (also referred to as microsatellites) can 
be found in the literature (Huang et al. 1999; Latch et 
al. 2002). Therefore we amplified and sequenced two 
turkey specific loci based on a study by Latch et al. 


(2002). While additional samples would be warranted 
to determine their utility, preliminary results suggest 
the WT75 locus varies within and between populations, 
and that the WT38-2 found little to no variation within 
and between populations. Small sample sizes confound 
interpretation of these results, yet may be of importance 
in future studies of Rio Grande wild turkeys in the High 
Plains (Table 1). 


Table 1. A comparison of the WT75 and WT38-2 loci microsatellites used in this study and Latch et al. (2002). 


Source Material 

Sample Identification 

WT75 locus Motif 

WT38-2 locus Motif 


TK136464 

[TGI + [TCG] 

[CA] 15 

Blood 


TK 136466 

[TC]„ 

[CA] 1S 


M3 07-2 

[GT] 5 + [CT] |0 

[CA] 15 

Egg Shell 


M307-4 

[GT]„ + [CT] 7 

[CA] 15 

Eastern wild turkey blood 

Latch et al. 2002 

[TC] lt 

[CA] 16 


Eggshell remains appear to be pragmatic sources 
of quality genetic material from wild turkeys. Despite 
low sample sizes, our data suggest no practical differ¬ 
ences between the quality of DNA isolated from blood 
and the quality of DNA isolated from egg shell remains. 


Further, the quality of DNA isolated from the chorioal¬ 
lantoic membrane of egg shell remains can be used in 
sex identification and other molecular techniques such 
as microsatellite amplification. 
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Appendix I 


List of specimens examined. Specimens were collectedfrom natural wild populations and are archived in the Natural 
Science Research Laboratory (NSRL) at the Museum of Texas Tech University (TTU). 


TK# 

DNA 

Concentration 

Sex 

Days from collection to 
suspension in lysis buffer 

140546 

NA 

Male 

73 

140547 

NA 

Male 

73 

140548 

NA 

Male 

73 

140549 

NA 

Female 

73 

140550 

NA 

Male 

73 

140551 

NA 

Male 

73 

140552 

NA 

Male 

73 

140553 

NA 

Male 

73 

140554 

NA 

Female 

73 

140555 

NA 

Male 

73 

140556 

NA 

Female 

73 

140557 

4.9 

Female 

107 

140558 

18.8 

Unknown 

107 

140559 

18.7 

Male 

107 

140560 

20.4 

Male 

107 

140561 

18.3 

Female 

107 

140562 

11 

Male 

107 

140563 

17.6 

Female 

107 

140564 

18.3 

Male 

107 

140565 

10.7 

Male 

107 

140566 

NA 

Male 

78 

140567 

NA 

Female 

78 

140568 

NA 

Male 

78 

140569 

NA 

Female 

78 

140570 

NA 

Male 

78 

140571 

NA 

Female 

78 

140572 

NA 

Male 

78 

140573 

NA 

Female 

78 

140574 

NA 

Female 

78 

140575 

NA 

Female 

78 

140576 

NA 

Male 

78 

140577 

66 

Female 

10 

140578 

170.1 

Female 

10 

140579 

166.6 

Male 

10 

140580 

117.7 

Male 

10 

140581 

125.4 

Female 

10 

140582 

99.6 

Female 

10 

140583 

126 

Female 

10 

140584 

178.8 

Male 

10 

140585 

90.7 

Female 

23 

140586 

39.2 

Female 

23 

140587 

74.2 

Male 

23 

140588 

139.9 

Male 

23 

140589 

121.3 

Female 

23 

140590 

145.8 

Female 

23 

140591 

163.6 

Female 

23 

140592 

134 

Male 

23 

140593 

169.1 

Unknown 

23 
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Appendix 1 (cont.) 


TK# 

DNA 

Concentration 

Sex 

Days from collection to 
suspension in lysis buffer 

140594 

227.2 

Male 

68 

140595 

63.9 

Male 

68 

140596 

108.5 

Male 

68 

140597 

335.6 

Female 

68 

140598 

77.6 

Female 

68 

140599 

111.2 

Male 

68 

140600 

161.2 

Male 

68 

140601 

437.2 

Male 

68 

140602 

302.7 

Female 

68 

140603 

33.3 

Male 

78 

140604 

27.5 

Female 

78 

140605 

44.7 

Female 

78 

140606 

45.9 

Male 

78 

140607 

26.8 

Male 

78 

140608 

19.4 

Female 

78 

140609 

19.8 

Female 

78 

140610 

23.9 

Male 

78 





